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Abstract. In this paper we study the motion of a charged particle on a 
Riemmanian surface under the influence of a positive magnetic field B. Using 
Moser's Twist Theorem and ideas from classical pcrtubation theory wc find 
sufficient conditions to perpetually trap the motion of a particle with a suffi- 
cient large charge in a neighborhood of a level set of the magnetic field. The 
conditions on the level set of the magnetic field that guarantee the trapping 
are local and hold near all non-degenerate critical local minima or maxima 
of B . Using sympletic reduction we apply the results of our work to certain 
5-'^-invariant magnetic fields on R^. 



1. Introduction 

The motion of a charge e on a Riemannian surface M can be described as 
a Hamiltonian dynamical system on T*M with metric Hamiltonian and twisted 
symplectic 2-form — d\ ei:* (3. Here A is the standard LiouviUe 1-form, 
TT : T*M M is the canonical projection, e is the particle's charge and /3 is a 
2-form on AI representing the magnetic field. A method for finding periodic orbits 



of this motion on a fixed energy level was introduced by V.I. Arnold (|Ar2, Ar3|) 



and developed by V.Ginzburg (|Gil, Gi2, Gi3]). The method relies on symplec 



tic topology techniques (Conley-Zehnder theorem |CZ|)or symplectic capacities) 
to prove the existence of some minimum number of periodic orbits. In this work 
we will combine ideas from classical pertubation theory and appropriately chosen 



changes of varibles (cf. [Ar| and |Gi2|) to study motion when the charge is large. 
The approximation involved with large charge is known as the guiding center ap- 
proximation. The motion of the particle is described as a fast rotation around a 
point that drifts slowly along a level line of the magnetic field. The point of rotation 
is known as the guiding center. Its behaviour is important since one is usually in- 
terested in confining charged particles to a bounded region of the phase space. We 
will find local conditions on a level set of a magnetic field on a surface with general 
Riemannian metric to perpetually trap the motion of the particle near this level 
set. By considering the twisted symplectic 2-form as the main element of our treat- 
ment we will write the system in a form suitable for applying Moser's Twist Map 
Theorem after applying successive diffeomorphisms. The motion will be trapped 
between the invariant circles predicted by the theorem. The non-degeneracy condi- 
tion of the theorem will be studied and neighborhoods where the condition of the 
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theorem holds will be characterized. Throughout the work the energy level will be 
fixed at E = ^. We state our main results: 

Theorem 1.1. Let M be a Riemannian surface. Let he the area form on M . Let 
w = dX + eTT*{Bfl), where A is the canonical 1-form on T* M , tt* : T* M M is 
the natural projection, B is a positive function on M and e > is a parameter. 
Denote by L^ = {B = c} a non-critical level set of the function B with the topology 
of an imbedded circle. Set j3 = Bfl = dA. Let F — |^Jj^ where Vi? denotes the 
gradient of B. If the level set Lc is such that 

(1) j iFd^ipdA^-^A^ 

then for any neighborhood N of Lc, there exists a number e* > such that for any 
e > e* the trajectories of the magnetic field problem with energy E = ^ and charge 
e will be contained on N for all times. 



Corollary 1.2. Let p G Al he a non-degenerate maximum or minimum point of the 
function B. Then there exists arbitrarily small neighborhoods N of p, and e* > 0, 
such that for all e > e* the integral curves of the magnetic problem with charge e, 
energy E — ^ and initial conditions in N will remain in N for all times . 



Corollary 1.3. Let s G M be a maximum or minimum submanifold of the function 
B with the topology of an imbedded circle. Then there exists arbitrarily small 
neighborhoods N of s, and e* > 0, such that for all e > e* the integral curves of 
the magnetic problem with charge e, energy E — ^ and initial conditions in N will 
remain in N for all times . 

The proof of the theorem will be based on a change of coordinates that will put 
the twisted symplectic 2-form edX + Bfl in the form 7 — e^dH^ A d9, where is an 
analytic function depending on e such that its first zero order term is given by -B^^ 
i.e, we will construct a diffeomorphism A such that A* {edX + B^) — j — e^dH^ A d9. 

The characteristic line field of this form is spanned by e^Xn, + where Xh^ 
is the hamiltonian vector field of iJe with respect to the area form 7. Thinking 
of the variable 9 as time allows us to introduce action-angle variables /, -0 on a 
neighborhood of the level set Lc of B, and to reduce the system to a time 27r map 
9 ^ 9 + 2Tr from this neighborhood to itself. The particular form of LL^ and 7 will 
allow us to invoke the twist theorem provided condition]^ holds for B. 

The paper is organized as follows. In section 2 we introduce the tools we need 
and describe the limit of large charge as a pertubative limit. In section 3 we 



construct the diffeomorphism A(theorem 3.1). In section 4 we apply Moser's twist 



theorem to our syst em and characterize the non-degeneracy condition. In section 



5 we prove theorem IT and its corollaries. In section 6 we apply our main theorem 
to a family of symmetric magnetic fields and finally in section 7 we prove some 
technical lemmas. 
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2. The Magnetic Problem for Large Charge 

Let M be a 2-dimensional oriented Riemannian manifold. Let /? = Bft be the 
magnetic field 2-form, where O is the area form on M and S is a non- vanishing 
function on M. The motion of a charge on the Riemannian surface M can be 
identified with a Hamiltonian flow (see [ pil[ |). Identifying the tangent and cotan- 
gent bundle by means of the Riemannian metric we denote by H : TM R the 
energy function H{X) = ^ where X G TM. Consider the Hamiltonian flow of 
H on T*M with respect to the twisted symplectic structure We = dX + en* 13. Here 
A = pidqi + P2dq2 is the canonical 1-form on T*M , and thus dX is the standard 
symplectic 2-form, tt* : TM ~ T*M M is the natural projection and e is the 
charge of the particle. Henceforth P and 7r*/3 will be identified. 

The twisted symplectic form We — dX + e(3 is not defined in the limit in which 
the charge is infinite. We will rescale the twisted 2-form so as to have an analytic 
limit. Let e — -. Set 

e 

= —We — — {dX + e/3) — edX + [3. 
e e 

Denoting by Xh the hamiltonian vector field of the magnetic problem we have 
dH = WeiXn, •) = ew^{XH, ■) = w^{eXH, •)• 

It follows that under the scaling Xh gets mapped to \Xh- This can be interpreted 
as a reparametrization of time. Thus if q{t) is a solution for X^, then q{^) is a 
solution curve for ^Xh- Therefore the scaling changes periodic orbits of period T for 
Xh to periodic orbits of period eT for \Xh- We thus have two associated problems. 
Given H we can consider the Hamiltonian problem given by the symplectic 2- 
form Wf, whose vector field we denoted by Xh or we can consider the Hamiltonian 
problem given by the symplectic 2-form w^^ whose vector field we denote by X^. 
From now on we will work only with Wf_. Write = \Xh- For future reference 
we collect these facts on the table 1.1. 



e = 1/e 


Expression 


Vector Field 


Time 


We 


dX + e(3 


Xh 


T 


We 


edX + P 


Xe = Xn/e 


eT 



TABLE 1.1 Relations between the scaled and non-scaled problem. 

The limit in which e goes to infinity can be understood as a pertubative limit. A 
closer look at Hamilton's equations will clarify this point. First we introduce some 
concepts necessary for our study. 

Denote hy S^M C TM the circle bundle over M with respect to the given 
Riemannian metric. Then S^M = {H = 5}- Let v be the geodesic vector field 
of M and v± be the vector field on S^M perpendicular to v (with respect to the 
orientation of the manifold). Denote by (f> the unit vector field in the direction of 
the fibers of the bundle ir : S^M ^ M with direction defined by the orientation 
of M. Thus, V, v± and form a global moving frame for S^M. Denote by v* , 
and 4>* the dual frame. Since energy is preserved the magnetic Hamiltonian vector 
field Xh at energy E = ^ is a vector field over S^M. The following lemma gives 
the structure equations for the circle bundle. 
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Lemma 2.1. The Lie-Brackets ofv,v± and (j) are given by 
where K is the Gaussian curvature of M . 



Proof. See |<N|. □ 
Dual to this lemma is 
Lemma 2.2. 

d(j)* = KvlAv*, dvl^v*A(j)*, dv* = (j)* Av*^. 
Remark: (p* is the connection 1-form for the Levi-Civita connection. 
Lemma 2.3. Let A denote the canonical 1-form on T*M. Then X = v* . 
Proof. In Darboux coordinates 



ij 



□ 



Corollary 2.4. Let J7 = w* A be the area form on M pulled back to S^M, A the 
canonical 1-form and v* , v*^_ and 0* as before. Then 

d{iy^n) = (j)* Av\^ dX. 
Consider the Hamiltonian equation for We, 

(2) w,{XH,-)^-dH. 
We write the Hamiltonian vector field as 

(3) Xh = v + ef^, 

i.e, as a geodesic component plus a fiber component. The function / is calculated 
by inserting this decomposition in equation ^, 

{dX + eBVL) {v + efcj)) = dH. 

Since , dX{v) = dH = (on SHd), n{v) ^ v\ and dX{<f>) = v\. We find 

that / = -B. Thus 

in S^M. Rescaling yields 



Xh 
Xh 



= ev 
e 



(4) 

Set X = X/eB and w — ^. Dividing by B we have that 

X = ew ~ (j). 

The integral curves of X^ and X are the same up to reparametrization. Thus the 
limit < e << 1 can be understood as the limit where the vertical vector field —(j) 
is perturbed by a 'small' vector field ew. 

This interpretation allow us to apply standard pertubation theory to eliminate, 
in the first approximation, the pertubation containing w by means of an appropriate 
change of coordinates. We indicate how to make this change of coordinates. Let 

xi — X{x) — ew{x) — 4>{x), 
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where / denotes derivation with respect to the new time. The basic idea is to puU 
back the vector field xi by the e-time flow $c of some vector field Y . For a small e 
we have (formally); 



s=0 



{xi) 



T = 



where t 



But 



d 



dr2 



^*x/ = [Y,xf]; 



$>/ = [Y, [Y, xr] 



and so on. Since x/ = ew — (j) the first order term of ^* (xf) is given by (recall that 



w 



The goal is to choose Y such that the first order term of the expansion is zero, i.e. 



(5) 



B 



- [Y, - 0. 



In classical pertubation theory, equation || is called the eikonal equation (see [ Ar|) 
Lemma 2.5. The vector field v±/ B solves the eikonal equation^. 
Proof. This follows directly from lemma 



2.1 



V .Wl,, V 1., - V V 

—,6] = — + — k>,v±] = = 0. 

B ^ B'^^ B B^^' ^ B B 



□ 



3. The Pertubative Approach 

The flow of the vector field v±/B gives the diffcomorphism that eliminates the 
first order term in the e expansion. In this section we prove the following theorem: 

Theorem 3.1. Assume Lc = {B = is a simple closed curve. For t sufficiently 
small (e > Oj, there exists a neighborhood Nc o/tt^^ (Lc) C S^M (depending on e) 
and a difjeomorphism A : Nc ^ Nc such that 

A*we = 7 - {e^dB-^ + e^dF) A de, 

where 7 is the area form on M , F is a function on S^M and 9 is a circular fiber 
coordinate on S^M . 

The change of coordinates A is involved. The basic idea is to construct it first 
on a section 9 = constant of the circle bundle and then extend it to the full circle 
bundle. The diffcomorphism on this section is built by pulling back the twisted 2- 
form by the fiow of the vector field v±/B introduced in section |[ In the flat case the 
triviality of the bundle implies that this is already sufficient to allow the application 
of the twist theorem. In the non-flat case we will need one diffcomorphism for each 
such section of the principal bundle. This apparent difficulty will be overcome by 
'gluing' together all these diffeomorphims. The main tool used to make this gluing 
will be a version of Moser's homotopy argument which says that two symplectic 
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2-forms that agree on a compact submaiiifold must agree on a neighborhood of the 
submanifold. The resulting diffeomorphism, now defined on the full bundle, will 
put the twisted symplectic 2-form in the final desired form. 

Proof. We start to build the diffeomorphism. We compute the pull-back of the 
twisted symplectic 2-form by the flow of v±/B. The choice of this vector field is 
dictated by its role in the pertubation problem as shown on the previous section. 

Lemma 3.2. The pull-back of Wc by ip^, the e -time flow ofv±/B is given by 

(6) ^,*w, ^Bn + ^^n- jdB-' A [de + /) + 0(e3), 

where 9 is the fiber coordinate function and f is a fiber-independent 1-form. 

Proof. This is a straightforward computation using Cartan's Formula and the struc- 
ture equations. The identity 0* = dd + f follows from the fact that (/)* is the 
Levi-Civita connection restricted to S^M. □ 



Corollary 3.3. The Flat Case : Suppose M is flat. Then K = 0, the bundle is 
trivial and we can take = dO i.e. / = 0. The pull-back computations reduce to 

Ljj± /gWe — —dX — edB^^ A dO. 
Ll±/BWe = {dB-^ ^eL,±/BdB-^) Ade. 
Ll^/B^e = [L,^/BdB-^ -eLl^/sdB-') Ad9. 

We obtain 

ipe*w^ = td\ + Bn + e{-d\-tdB-'^ Ade) + — {dB-^ -eL„±/BdB-^) Ade + 

(7) i- (L,^/BdB-^ ~ eLl^/gdB-^) Ad0 + ... 

which simplifies to 

2 3 4 

ip,*We ^Bn- -^dB-^ Ad9- —L„^/BdB-^ A d9 - —Ll^j^dB'^ Ade + ... 
giving 

ip,*w^ ^Bn- e^dH^ A d9, 

where 

?— 1 ^ ^ 

So, for the planar case, the first diffeomorphism is already sufficient to bring the 
twisted symplectic 2-form to the desired form. One should compared this expression 
with the one found by Littlejohn Q. 

A second diffeomorphism will be necessary if ^ 0. Denote We = (pe*We. First we 
will restrict to the level set 9 = constant. In this case we have that d9 = and 
(||) becomes 

(8) w,^Bn + ^^n-^dB-^ Af-h... 

Z 13 Z 
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Lemma 3.4. For a constant 9 and a sufficiently small e (e > 0) there exist neigh- 
borhoods No and Ni (depending on e) of Tr^^{Lc) C S^M, and a diffeomorphism 
^e,e : No ^ Ni such that ^* gW^ — BVL. Moreover 

- /rf + e^C(e) 
with C(e) uniformly hounded in the norm. 

Proof. 7r~^(Lc) H {0 = constant^ has dimension 1, consequently the forms We and 
BVt degenerate. Since Lc is compact it follows from theorem of Weinstein (a version 
of Moser's homotopy argument) that there exists a diffeomorphism $£.9 such that 
gW^ = Bfl. The order of the diffeomorphism is given by tracing back Moser's 
homotopy argument. This is done in section ^ □ 

Composing the two diffeomorphims ^^,9 and ipe we obtain on the section 6 = 
constant of the principal bundle that 

For a fixed e, ^e.e is a family of diffeomorphisms parametrized by 0, each defined 
only on its own constant section of the principal bundle. We will paste all those 
diffeomorphisms together. Let p G S^M . With respect to our local trivialization we 
write p = (m, 9) for an arbitrary point of S^M in the given neighborhood. Define 
: S^M S^M by 

S,(m,0) = ($,,e(m),0). 

Lemma 3.5. Assume = {B — is a simple closed curve. For e sufficiently 
small (e>Q), there exists a neighborhood Nc of Lc (depending on e) and a diffeo- 
morphism S : Nc Nc such that 

E*w, = BQ- {e^dB-^ + e^R) A d9, 

for some 1-form R. 

Proof. Since restricted to 9 — constant is equal to ^^.e if follows that S^Wg 
differs from ^* gWe by a factor of the form a A d9, for some 1-form a, i.e., we can 
write 

(9) E*We^ Bn + aAd9. 

The one-form a is defined mod d9 and can be taken to be 



For any diffeomorphism S we have that i5* {ixct) = is,x iS*a) (see e.g. |AM 
Therefore 



(10) = 

_d_ 

' 00 



So we must calculate 5e^^. We do this on the following lemma. 



Lemma 3.6. 

— = 1- e z 

" 89 39 

where z is a horizontal vector field relative to our trivialization. 

Proof. The proof of this lemma is given in section |^ □ 
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Using this lemma we compute ( p^ ) , 



(11) i a E*We — S*(z- a We) = We) + — e*(izWe)- 

We recall the expression for (formula ^ 

£2 ^ 2 

zDe = ^*We = Bn + ——n - —dB-^ A {d9 + /) + 0(6^). 
2 ij 2 



It follows that 



iaw, = -—dB-^ - e^G, 



for some 1-form G. We have that ^ can be writen as 



Using lemma 3.4 we have 



i^we = iU + e'G*ie, 6)){--dB-^ + e^G) + e^S*(i,tSe), 



givmg 



i a Wt ^ -—dB ^ - e^R, 
for some 1-form R. Finally this implies 

(12) E^We = Bn + aAde = BQ- {e^dB-^ + e^R) A dO 



which proves lemma 3.5. □ 



The 2-form (12) is almost in the final form. The 1-form R may not be necessarily 
exact. A final change of coordinates is needed. So far we have been working on a 
neighborhood of a constant level curve {B — c} of the magnetic field. By hypothesis 
this has the topology of an annulus. We introduce Fermi coordinates {x, y) on this 



annulus (see e.g. |Mc]). The Fermi coordinates parametrize the Lc neighborhood 
A^o on the following way: Given a point m G Nq, consider the segment of geodesic 
S that connectes m to Lc and has minimal length. Let y denote its oriented 
length. Denote by x the oriented length from some arbitraty point xq in Lc to the 
intersection of S with Lc- Using these coordinates we define 



aie) = / Rdx, 
Jb=c 



where l{c) denotes the length of the level set Lc- 

Lemma 3.7. With a{9) as above there exist functions F and g such that 

R = dF + a{d)dx + gdO. 

Proof. Since B is 9 independent and is closed (|l^) implies that 

(13) dRAde = 0. 

Let 

R = R^dx + Rydy + Rgd9, 
and denote by R* the restriction of i? to a constant section of S^M, i.e. 

R* = i?|e=c = Rx{x,y,c)dx -f Ry{x,y,c)dy. 
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Thus (|13D implies that 



f dRy dRx 



dRhde^ ~ ] dx Ady Ade = 0, 

\ ox ox J 

from which we see that R* is closed. Now define a 1-form ^ on the section 6 = chy 

(14) ^ = i?* - a{c)dx. 
Thus 

(15) = and / C = 0, 
giving that 

a(c) ^ I R* 



for any L' homotopic to Lc (since R* is closed). Since the section {6 = constant} 
has the homotopy type of an annulus we have that equations |l^ together imply 
that 

for some function F. Thus we can write |l^ as 

R* =dF + a{c)dx. 

But R* difers from i? by a factor of the form gdd for some function g and we can 
write that 

R = dF + a{e)dx + gdO, 
proving the lemma. □ 



Using lemma 3.7 we can rewrite [12| obtaining that 

(16) S*w;e =Bn- (e^dB-^ + e^dF + e^a{0)dx) A dO. 



Define a new variable 



and the diffeomorphism 



y ^ I Bdy. 



X^x, 



(17) \ Y^y-e^J^a{e)de, 

e = 61. 



Substituting back in 16, and denoting by the pull-back of by S* we obtain 
that 

= dX A dY + e^a{Q)dX AdO- {e^dR-^ + C'dF + e'a{Q)dX) A dQ, 

i.e, 

Mv^ ^ dX A dY - {e^dR-^ + e^dF) A dO, 

that we write as 

(18) w.^dXAdY + e'^dH.AdQ, 

where = B^^ + eF which proves theorem 3.1 □ 
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4. Moser's Twist Theorem 



The characteristic hne field of (|18|) can then be thought as being generated by the 
Hamiltonian system given by the Hamiltonian H,, and symplectic 2-form dX A dY. 
We introduce action angle variables /, ■0 on a neighborhood of the level set Lc of 
B. Using the Fermi coordinates (see section^ we have that the level set is given by 
Lc — {y = 0,9 = 0} and is contained in an open set TV = {—S < y < 5,9 — 0} for 
some S > and N C Nq. Denote by N the image of under the diffeomorphism. 
Thus we have that N becomes N={--|<r<-|,e = 0}. Define 

(19) I{X,Y) ^ I dX AdY. 

where A{X,Y) is the annulus y<F<OorO<F<y. Calling ip the variable 
canonicaly conjugate to / we have that 

H,= '^B-\l) + Oie'), 

dX AdY ^dl A d%l). 

We consider the Hamiltonian equations with on A{X,Y). By integrating and 
taking the time {9 — 2tt) flow we get 



(20) 



/i =/ + 0(e3), 

! d 

At this point we recall the statement of Moser's Twist map theorem: 



V'l -^ + ^e2^S-i(/) + 0(e'). 



Theorem 4.1. (Moser, 1962 Mos|J Let ^ be a measure preserving map <&(i?, 9) 

{Ri, 9i) given by 



(21) 



Ri^R + e^f{R,9,e), 



71 



9 + a + e''-f{R) + e^g{R,9,e), 



where "f/(R) 7^ 0, k < I, f and g bounded. Then for e sufficiently small there exists 
an invariant curve T surrounding R — 1. More precisely there exists a differentiate 
closed curve 

^ ^ \ ^? = </) + G(</.,e) 

with F, G of period 2t: in (j) which is invariant under the map 



To apply the twist map theorem 4.1 to ^it is sufficient that -^B~^{I) 7^ at 
Lc. This is the so called non-degeneracy condition and plays an essential role in 
our problem. In what follows we will characterize the non-degeneracy condition. 

4.1. The non-degeneracy condition. 

Definition 4.2. A level set Lc — {B = c} is called non-degenerate if 
along Lc. 
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Using the chain rule and assuming is a non-critical level set of B we can write 
that ( psj ) is equivalent to 

(24) = -i^. 

^ ' dB^ BdB 

Form ( |2^ ) is particularly useful since one can explicitly calculate the derivative of / 
with respect to B. To do so we first note that /, the action variable, can be easily 
expressed in terms of the Fermi coordinates x and y. Since 



I{X, Y)= ( dX AdY 

J A{X,Y) 



we have that 



/(x, y) — Bdx A dy. 

J A(x.v) 



I A(x,y) 

where A{x^y) is the preimage of A{X,Y) under the diffeomorphism (|1 
Lemma 4.3. 

where ds is the arc length of Lc- 

Proof. Choose a vector potential A = A^dx + Aydy for the 2-form Bdx A dy. Thus 
we have by Stokes theorem that 



I{B) = I A. 



Now consider the vector field F — \vb\'^ • ^'^^ time t flow. Thus 

implying that (p^B ^ B + e. We can write that 



I{B + e)^ f A^ I 



\A. 



Differentiating with respect to e and using Cartan's formula we obtain that 

d_ 



-I{B + e)= {dipA + iFdA). 



The first term on the integration vanishes since dip A integrates to zero along the 
closed curve L^. To get the other term, observe that for any smooth simple closed 
curve with smooth normal n we have indxAdy = ds the arc length. Since F = j^^^n 

and dA — Bdx A dy we find that ipdA — j^^ds as desired. □ 
Corollary 4.4. 

(26) dB^ = J ^pdiipdA). 



Proof. Since the second derivative is expressed as the integral of a 1-form, the same 
reasoning applied to compute the first derivative can be applied to compute the 
second. □ 
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The non-degeneracy condition can be written in a more explicit way by noting that 
ds can be written as f{x, y)dx for some function / > 0. A long and straightforward 
computation give us that the level Lc is nondegenerate if 

^''^ 1 - ^ < '''' ^) ^ W < 4 ^ 

This expression for the non-degeneracy condition will be used to prove corollary 



1.2 



5. Proof of Theorem [l] 

Proof. We look to the magnetic problem given by the 2-form = ed\ + BVl for 
H — ^2 where e = i. According to table 1.1 the Hamiltonian vector field for is 
given by a scaling of the Hamiltonian vector field for w. Thus their characteristics 



differ only by a time reparametrization. By theorem 3.1 there is a neighborhood N 
of Lc and a diffeomorphism : N ^ N such that 

S>, = 7 - (e^dS-i + e^dF) A dO. 

Thus the characteristic line bundle of S*ii'e is spanned by the Hamiltonian vector 
field given by the Hamiltonian system with 

= e^B-^ + e'F 

and symplectic 2-form Vi (where by abuse of notation we denote S*ri by fi). In 
section ^ we introduced action- angle coordinates on the neighborhood of "B* N and 
reduced the dynamics of the Hamiltonian system (ffe, fi) on this neighborhood to 
the dynamics of a twist map. The nondegeneracy condition of Moser's twist theorem 
is calculated to be equation |l| according to ( p^ ) and lemma 4.3. So Moser's twist 



theorem applies and invariant circles exist. By dimensionality, they trap the charge 
for all times t > □ 



Proof. ( of Corollary [LJJ Observe that each term in (|27| ) depends on diferent powers 
of \VB\. A careful analysis of this fact (done in section |^) implies that ( p7| ) can 
not be satisfied near a critical point of B. So for a sufficiently large e we can find 
a small neighborhood of p where condition (E3[) is satisfied everywhere. □ 



Proof, (of Corollary l.c) The proof is by contradiction. Suppose that given a 
neighborhood TV of Lc we have that all the level sets contained on N are degenerate. 
This implies that exists a constant K such that 



on N. This is equivalent to 
What implies that 



1 dB -1 d 
WdT ^ 'B^'dB 



K. 

dl r B 



ds. 



KB^ dB J^^ |VB| 
Since the right hand side is positive, the constant K must be negative and writing 
\K\ = it follows that 

\K\ f I 1(c) 
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where l{c) denotes the length of the level Lc- Thus we have that 
(28) HO < 1^1— r^l). 



where m is the minimum of B on N, the closure of N. (|28|) implies that l{c) 
is dominated by max {\VB(c)\). Since N was arbitrary this implies that 1(c) the 
length of the critical level set must be zero, but this contradicts the hyphothesis. □ 

6. The Symmetric Case 

In this section we deal with symmetric magnetic fields on a three dimensional 
manifold, namely M.^ with the z-axis removed. Using symplectic reduction the 
dimension of the system will be reduced allowing us to apply the main theorem of 
this work. The non-degeneracy condition of one of the components of the magnetic 
field will be sufficient to apply Moser's Twist theorem. 

Let M be the standard euclidean three dimensional space with the z-axis re- 
moved. Let act on R"^ by rotation around the z-axis. Let an invariant 1-form 
A on M be given. Using cylindrical coordinates we write 

A = Ardr + Agde + A^dz. 

The magnetic system given by /? = dA and the purely metric hamiltonian is equiv- 
alent to the Hamiltonian system given by 

(pe-eAg)^ {Pr-eArf {p, - eA,f 
+ ^ + ^ 

and symplectic 2-form 

w — dpe Ad9 + dpr Adr + dpz A dz. 
The momentum map for rotation around the z-axis is pg . Fixing pg — AI it follows 



by a theorem of M. Kummer |Ku| that the reduced Hamiltonian is 



where 

[M eAef 
V<^ff - ^2 . 

the reduced symplectic 2-form is given by 

/ dA dA \ 

U — dpr A dr + dpz A dz + e [ ' — ] dr A dz 

\ or oz J 

and the reduced manifold is given by = M x M"*" with coordinates (z, r) which 
can be visualized as a half-plane with boundary the z-axis. Writing (3 as 

P = Bgdr Adz + Brdd Adz + B^dr A dO 

we see that 

= dpr A dr + dpz A dz + eBgdr A dz = dX + eBgdr A dz 
here dA is the canonical 1-form on T*N. Assuming that 

E > Veff 

we have that the hypersurfaces 

Me = [h = E^ 
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and 



M,^{(£-K.„)-(|4).l} 



are equal. It follows that the closed characteristics of f^|jv/E ^'^^ equivalent to the 
closed characteristics of ^\mi- Thus we have that the reduced dynamics is (up to 
reparametrization) given by the purely metric Hamiltonian 



H = 



2 

'as 



2 



(where p = {pz,Pr) and (?_e is the metric gE — {E — Veff)~^g where g denotes the 
Euclidean metric on T*N) and symplectic 2-form fl. Assume that a level set of 
B = (Bg, Bz, Bf) restricted to is a simple closed curve, call it Lc, and that for 
a neighborhood of Lc, the non-degeneracy condition holds, that is to say that Bg 
satisfies (|7|) on a neighborhood of Lc- Then we can apply Moser's twist theorem 
as in the proof of our main theorem. 



This result should be compared with the works of F. True and M. Braun |Br|. 
They considered symmetric magnetic fields on M.^ given by vector potentials of the 
type Agdr A dz, which is to say Bg = 0. Those are extremely degenerate magnetic 
fields in the sense of (|2^) and so our work does not apply. Also in their case, 
the absence ot the Bg component of their magnetic fields simplified the problem, 
after the reduction of the symmetry, to a problem of the type kinetic plus potential. 
Braun's work is, to the author's knowledge, the first to apply Moser's twist theorem 
to a magnetic problem. Truc's work is remarkable since she was able to prove the 
trapping of the particle in the case where the magnetic moment was not convex. To 
our knowledge we are the first to generalize the problem to a Riemmanian surface 
and to deal with problems where the magnetic field could not be eliminated by the 
reduction of the symmetry. It is in this sense that our work complement theirs 
since we dealt with disjoint classes of magnetic field. 



7. Technical Proofs 



Proof, (of Corollary [3j|j We follow McDuff fMcDj]. The first part of |3J is an 
imediate consequence of Moser's homotopy argument. The diff'eomorphism $j g is 
built by realizing it as the flow of a vector field X on M. First we consider the 
family of 2-forms Wt ~ (1 — t)Bil + tw^ for t E [0, 1]. We want a vector field X such 
that its flow <I>e,e is such that 

^Ig^.wt = Bn. 

(in what follows we will omit the 9 and e dependency from the notation). Differ- 
entiating this relation we obtain that 

+ = 0, 

i.e. that, 

$* {Lxwt + [we - Bn)) = 0. 
Since $ is a diffeomorphism and Wt is closed it follows that this can be satisfied if 
(29) d{ixwt) + {w, ~ Bn) = 0. 
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We want to solve equation for X. Since Wt is non-degenerate, it suffices to find 
a family of 1-forms at such that 

(30) 0-4 1^,^^=0, d(Ji=w,-Bn. 
In fact, conditions (^ ) imply 

(31) CTt + ixt^f = 0. 

To construct at we consider the restriction of the exponential map to the normal 
bundle TLjr of the submanifold Lc with respect to the Riemannian metric on M . 
We denote this restriction by 

(32) exp : TL^ M. 
Consider the neighborhood of the zero section 

J7e = {{m,v) e TM I TO e M, -y G TL^, \v\ < e} . 

Then the restriction of the exponential map to Ue is a diffeomorphism onto Nq = 
exp{Ue) for e > sufficiently small. Let {p,v) £ TM. Define 7r2 as 

TT2{p,v) = V. 

Now define ipt : Nq ^ Nq for < t < I by 

^til) = Gxp{q,tv), where v — 7:2 (^exp^^ ((7)) . 

Then ^t is a diffeomorphism for t > and we have ipoiNo) C Lc, ■01 = ^d,, and 
V'tUc = Calling T = We — BCl we define 

(33) --/l^**- 

It's a direct computation to check that the family of 1-forms at has the desired 
properties. Thus the family {at) satisfy conditions ^ and it follows that ^e,e the 
flow of X, the vector field that satisfies ( pl| ) is such that 

^Igw^ = BQ. 

It follows that 

where R stands for the tail of the expansion of r. Denoting F = (^f^ — — ) we 
have 

that we write as. 



dt jr. dt 



at ~ e'^ao,t + e^ci^t 



where ctq,* = /q ■^ip*F and ci^t = /g ^ijfR. To solve equation ( p9| ) we write 



where X^^j is the solution of 



X — e^Xo^t + e'^^i.t, 
CTi,* -I- ix.^tW = 0. 
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for i = 0, 1. We write X = e^Y where Y = X^^t + ^Xi^f Since X is a continuous 
vector field it follows that |y| is bounded on a neighborhood of A^o by a constant 
K. The equation for the flow <i> is 

Chosing $(0) = po € Nq and integrating it follows 

$(po,t)-$(po,0)= / X{Mt))dt. 

Jo 

Writing <1> the flow of X as $ = /d + e^C(e) we obtain 

PQ + e'C (po, e, t)-po = e^ f Y(Mt))dt, 

Jo 

giving 



|e'C(e)| < £2 / \Y\dt < e^K 



and the result follows. 



□ 



Remark 7.1. An important consequence of this lemma is that Xo^t is fiber indepen- 
dent. This follows directly from the definition. 



Proof, (of Corollary S.t) Using the Fermi coordinates let p = (x, y, 6) be a point of 
S^M . Omiting the e dependence of 2 for notational convenience we can write that 

S(p) = (S.(p),S,(p),0). 

The derivative of S can be calculated as 

/ 35^ d'Eso c?5x \ 
< dx dj/ 99 \ 

d^y d^y 

dx dy do 





which give us that 
(34) 



1 / 



1 



where p/ = ^(p). Recalling the definitions and lemma 3^ we can write that 



^2 dCrr 

,2 OCy 
09 ■ 



Since ^^^g = Id + e^C(e, 9) is the flow oi X = e'^Xo^t + £^Xi,t we have 



dt 



Integrating this equation we have that 

e2C(e,0)= / {e^Xo,i + e''X,,t)o{U + e^Cie,9))dt, 
Jo 

that simplifies to, 

(35) C{e, 9)^ f {XoAU + e'C(e, 9)) + eX,,t + e'^i.t(C(e, 9))} dt. 
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Expanding Xo^ti^d + e^C(e, 6*)) in e, inserting the resulting expression in ( |35| ) and 
collecting terms we have 



C{e, 9)= f Xo^tdt + e f XrcH, 
Jo Jo 



for some X^. Since Xq t is 6 independent, we obtain 

This is to say that ^'^gg'^'* is of order at least one in e and we write that 



where 




1 



" / N ^ 9 3 



where 



proving the lemma. □ 




Proof, (of Corollary LJ j We will study each term of ( p7[ ) individually. Let B{p) — 
a. Since p is a nondegenerate minimum (or maximum) point we have by Morse's 
lemma that there is a neighborhood V and a system of coordinates {u, v) on V such 
that B{u, v) — (u^ + u^) where the plus sign is used if the point is a minimum 
and the minus sign is used if the point is a maximum. 

Lemma 7.2. Let r^ = + and let r^ = \ {u, Then we have that 
(36) |v5|2 = ^«^^[^->'-^' 



9 

The proof of this lemma is straightforward. Note that (pq) implies that trace[gi. 



r: 



2 



> on y — {p}. For a level set c not equal to a and such that Lc CV ~ {p} we 
have that rg is constant (since B is) and we can write that 

5|VB|2 > frace[g^j]rl, 

i.e., 

\ 9 

Now we can see that the first term of ( p7| ) is bounded, in fact we have that 

3 fdx < 3 (—^—-] ' / fdx, 



\yB\ \trace[g. 
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but on V we have that g and / are bounded above and trace[gij] is bounded below. 
So we can find constants M and K such that 

M > ( — r] \ and K>f 



\ trace [g^j j / 
on y, allowing us to write that 

(37) / T^fdx < ^ f dx^ ^27rre = 2ttMK. 

This give us that the first term of is bounded for any level set close enough to 
p as claimed. Now we look to the second term 

(38) j^<V\\/B\,^B>fdx. 

For a constant metric one computes that 

< V|VS|, VB >= |VB|. 

381) becomes 

^ :Jdx, 



considering lemma ( ^6[ ) , and proceeding as before we have that 
39 / T^T7^fdx>— -^dx>—^; 

The case of a nonconstant metric can be reduced to the case of a constant metric 
by choosing normal coordinates on a neighborhood of p. 
The last term of (p7|) can be minorated by noticing that 



< V/, VB > dx 



VB 



< I ^|V/|<G ^ ' 



for some constant G (since |V/| and B are bounded on V). Now observing that 



|VB|2 = r2, (^ implies 



And we can write 



1 H 
< — 



for some constant H such that 



onV. It follows that 



|VB|2 

g + 1 

trace[gi- 



< H 



(40) / — ^<V/,VB>rfx<27rG-. 

Jl,_ |vB|3 re 

Now considering (|3^), (||), (|^) we see that (|2^) implies 

c K H 
MK>——-G—, 

M rl Te 

which can not be satisfied if rg is small enough. Thus for a sufficiently small 
neighborhood of p all the level sets are non-degenerate. □ 
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